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Abstract 

In the framework of the symplectic extension of the Interacting 
Vector Boson Model (IVBM) a good description of the first excited 
positive and negative parity bands of the nuclei in the rare earth 
and the actinide region is achieved. The bands investigated in the 
model are extended to very high angular momenta as a result of 
their consideration as "yrast" bands with respect to the symplectic 
classification of the basis states. The analysis of the eigenvalues of 
the model Hamiltonian reveals the presence of an interaction between 
these bands. Due to this iter action the AL = 1 staggering effect 
between the energies of the states of two bands is also reproduced 
including the "beat" patterns. 



1 Introduction 

The existence of nuclei with stable deformed shapes was realized early in 
the history of nuclear physics. The observation of large quadrupole moments 
led to the suggestion that some nuclei might have spheroidal shapes, which 
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was confirmed by the observation of rotational band structures and mea- 
surements of their properties. For most deformed nuclei, a description as an 
axial- and reflection-symmetric spheroid is adequate to reproduce the band's 
spectroscopy Because such a shape is symmetric under space inversion, all 
members of the rotational band have the same parity. However, with the first 
observation of negative parity states near the ground state, the possibility 
arose that some nuclei might have an asymmetric shape under reflection. 

On the other hand, whenever symmetry breaking appears new behavior 
of the many-body system is expected. Reflection symmetry breaking is as- 
sociated with a static octupole deformation which is expected to determine 
new collective features for the nuclear system. 

Extensive investigations into the structure of nuclei with low-lying neg- 
ative parity states has led to the conclusion that, while reflection asym- 
metric shapes can play a role in the band structure, they are not as stable 
as the familiar quadrupole deformations. The rotational spectra of some 
even-even nuclei in the rare earth and light actinide region exhibit, next to 
the ground band, a negative parity band which consists of the states with 
I* = l~,3 - ,5~, ... . These two bands are displaced from each other, which 
means that fluctuations back to space symmetric shapes must also be signif- 
icant. Experimentally the presence of "octupole" bands for some isotopes 
from the light actinide and rare earth region pQ is firmly established. 

There is a large variety of models that try to describe this behaviour of 
the low-ling states of deformed nuclei. Particularly successful are algebraic 
models based on symmetry principles. The introduction of an additional 
octupole degrees of freedom is a common feature of most of those models . 

The prescription for describing negative parity states by the addition of 
an / boson to the usual s and d of the IBM was first mentioned by Iachello 
and Arima j2]. It was suggested |3| that the inclusion of a p boson to the 
s, d and / bosons may play an important role in the description of these 
collective states. 

The coherent state method (CSM) was applied by Alonso et al. to the spdf 
SU(3) Hamiltonian with quadrupole and octupole interaction Recently 
A. A. Raduta and D. Ionescu jHj have used a generalization of the CSM . 
They suggested that both ground and octupole bands may be considered 
as being projected from a single deformed intrinsic state that exhibits both 
quadrupole and octupole deformations. 

Another collective model based on point symmetry group considerations 
jU] has also been used very successfully for the description of the energy 
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levels of the ground and octupole bands and reproduces odd-even staggering 
between these levels j7]. In this model the octupole field is parametrized by 
irreducible representations of the octachedron point symmetry group. 

The introduction of an octupole degrees of freedom in the presence of 
comparatively large number of free parameters in all of these models allows 
for the reproduction of the experimental data of the energies of the negative 
parity states, at least in the low spin region. 

In the beginning of the 1980's a phenomenological algebraic model called 
the Interacting Vector Boson Model (IVBM) was introduced jH] • This model 
is a generalization of the phenomenological broken SU(3) symmetry model 

which provided a good description of the low- lying ground and 7 bands 
[TUj of well deformed even-even nuclei. Its advantages were incorporated 
into the rotational limit of the IVBM ^T] , with an a good description of all 
the positive parity bands of nuclei in the rare earth and the actinide region. 
Moreover, the £7(6) extension of the model contains such sequences of SU(3) 
multiplets, some of which proves to be convenient for the description of the 
low-lying negative parity bands [T2| . 

With the recent advance of the experimental technique the investigated 
bands were extended to very high angular momenta PQ. This motivated a 
new approach within the framework of the model aimed at a description of 
the first positive and negative bands, up to very high spins. In this new 
application, we make use of the symplectic extension of the IVBM |T3] . This 
allows these bands to be considered as yrast bands in the sense, that we take 
into account the states with a given L, which minimize the energy values 
with respect to N. N is the eigenvalue of the total number of bosons that 
build the basis states of the model. Its eigenvalue changes as AiV = 2 in the 
infinite spaces of the boson representation of Sp(12,R). When considering 
the dynamical symmetry of the symplectic extension of the model through 
the maximal compact subgroup £7(6) D Sp(12,R), we obtain the exactly 
solvable rotational limit with a Hamiltonian, diagonal in a basis defined by 
the irreducible representations of the corresponding chain of subgroups. The 
measured energies of the ground and octupole bands in even-even nuclei 
from the rare earth and actinide regions are reproduced in the model with 
rather good accuracy. The analysis of the obtained results shows that this is 
due to the appearance of a vibrational type term that influences the yrast 
energies. This term also plays the role of an interaction between the two 
considered bands, and is the reason for the correct reproduction of the 
odd-even staggering of their energies. 
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2 Algebraic Basis of the IVBM. 



We start with a brief review of the model's assumptions and definitions. The 
IVBM is based on the introduction of two kinds of vector bosons (called p and 
n bosons), that "built up" the collective excitations in the nuclear system. 
The creation operators of these bosons are assumed to be SO (3) vectors and 
they transform according to two independent fundamental representations 
(1,0) of the group 577(3) . These bosons form a "pseudospin" doublet of 
the U(2) group and differ in their "pseudospin" projection a = ±|. The 
introduction of this additional degree of freedom leads to the extension of 
the SU(3) symmetry to U(6). Then the operators 

u m( a =\)=PL u i( a = -\) = n i (!) 

transform according to the fundamental representation [l]g of the group 
U(Q). The annihilation operators are obtained by the conjugation (it+(a;))t = 
u m (a) and transform according to the conjugate SU(3) representations (0,1). 
The bilinear products of the creation and annihilation operators of the two 
vector bosons generate the noncompact symplectic group Sp(12,R) 

G L M (a,(3) =% m C^ m Ma)u m W), 

A L M (a,(3) =f m C<W%W> (2) 

where C^ m are the usual Clebsh-Gordon coefficients and L and M define 
the transformational properties of (J2J) under rotations. 

We consider 5^(12, R) to be the group of the dynamical symmetry of the 
model [S]. Hence the most general one and two-body Hamiltonian can be 
expressed in terms of its generators . Using commutation relations between 
the F^(a,(3) and G^(a,/3), the full range of number of bosons preserving 
Hamiltonian can be expressed in terms of operators A^(a,^): 

H = Y,ho(a,p)A°(a,l3) + £ E (-l) M V L (a(3; yS)A L M (a, i)A L _ M ((3, 5), 

a, /3 M,L af3~f8 

(3) 

where h (a,j3) and V L (aj3; j5) are phenomenological constants. 



4 



Being a noncompact group, the representations of Sp(12, R) are of infinite 
dimension, which makes it impossible to diagonalize the most general Hamil- 
tonian. The operators A\j{a,f3) generate the maximal compact subgroup of 
Sp(12, R), namely the group U(6): 

Sp(12,R) D U(Q) 

So the even and odd unitary irreducible representations (UIR) of Sp(12,R) 
split into an infinite but countable number of symmetric UIR of J7(6) of the 
type [N, 0,0, 0,0,0] = [N] 6 , where N = 0,2,4, ... for the even set (see Table 
1) and iV = 1,3,5,... for the odd set [T3j. These subspaces are of finite 
dimension, which simplifies the problem of diagonalization. Therefore the 
complete spectrum of the system can be calculated only trough the diagonal- 
ization of the Hamiltonian in the subspaces of all the UIR of U(6), belonging 
to a given UIR oiSp(12,R). 

The rotational limit jTTj of the model is further defined by the chain: 

U(6) D SU(3) ® U{2) D SO(3) ® U(l) (4) 

[N] (A, n) (N, T) K L T (5) 

where the labels below the subgroups are the quantum numbers (JHJ) corre- 
sponding to their irreducible representations. Their values are obtained by 
means of standard reduction rules and are given in ^T]. In this limit the 
operators of the physical observables are the angular momentum operator 

L M = -V2~Y1 A \A a ^) 

a 

and the truncated ( "Elliott" ) quadrupole operator 



Qm = Vq a m 

a 

which define the algebra of SU(3). 

The operators of the "pseudospin" components and the number of 
bosons N : 



T = -J-[A°(p,p)-A°(n,n)}; N = -V3[A°(p,p) + A°(n,n)], 
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define the algebra of U(2). 

Since the reduction from 17(6) to SO (3) is carried out by the mutually 
complementary groups 517(3) and 17(2), their quantum numbers are related 
in the following way: 

T = ± N = 2fi + X (6) 
Making use of the latter we can write the basis as 

| [N) 6 ; (A, n); K, L, M; T > =| (N, T); K, L, M; T ) (7) 
The ground state of the system is: 

| ) =| (N = 0, T = 0); K = 0, L = 0, M = 0; T = ) (8) 
which is the vacuum state for the 5p(12, R) group. 

2.1 The symplectic extension of IVBM 

The basis states associated with the even irreducible representation of the 
Sp(12, R) can be constructed by the application of powers of raising gener- 
ators F^(a,/3) of the same group. Each raising operator will increase the 
number of bosons iV by two. As a result we get a realization of the reduction 
scheme |14j : 

Sp(12,R) N>U(6) SU{2)xSU{3) SU{3) (9) 

The Sp(12,R) classification scheme for the £77(3) boson representations for 
even value of the number of bosons iV is shown on Table 1. Each row (fixed 
N) of the table corresponds to a given irreducible representation of the 17(6). 
Then the possible values for the pseudospin are T=y,y — 1,... and are 
given in the column next to the respective value of N. Thus when iV and T 
are fixed, 2T + 1 equivalent representations of the group 517(3) arise. Each 
of them is labelled by the eigenvalues of the operator T : — T, —T + 1, T, 
defining the columns of Table 1. The same SU(3) representations (A, y) arise 
for the positive and negative eigenvalues of T . 

Hence, in the framework of the discussed boson representation of the 
Sp(12, R) algebra all possible irreducible representations of the group 517(3) 
are determined uniquely through all possible sets of the eigenvalues of the 
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Hermitian operators N and T 2 . The equivalent use of the (A,//) labels fa- 
cilitates the final reduction to the SO (3) representations, which define the 
angular momentum L and its projection M. The multiplicity index K ap- 
pearing in this reduction is related to the projection of L in the body fixed 
frame and is used with the parity (tt) to label the different bands {K 71 ) in the 
energy spectra of the nuclei. We define the parity of the states as n = (— 1) T . 
This allows us to describe both positive and negative bands. 

2.2 The energy spectrum 

The Hamiltonian, corresponding to the considered, rotational limit of IVBM, 
is expressed in terms of the first and second order invariant operators of the 
different subgroups in the chain Pjl: 

H = aN + a 6 K e + a 3 K 3 + a x K x + (3 3 tt 3 , (10) 

where K n are the quadratic invariant operators of the U{n) - groups, 7r 3 is 
of the 50(3) second order Casimir operator. As a result of the connections 
(JHJ) the Casimir operator K 3 with eigenvalue (A 2 + fi 2 + A/i + 3A + 3/i), is 
expressed in terms of the operators iV and T: 

K 3 = 2Q 2 + -L 2 = -N 2 + N + T 2 

Making use of the above relation, Hamiltonian ()10|) takes the form 

# = aAr + ^ 2 + a 3 T 2 + /Vr3 + a 1 T 2 , (11) 

and is obviously diagonal in the basis (J7J) labelled by the quantum numbers of 
the subgroups of group-subgroup chain (j3J). Its eigenvalues are the energies 
of the basis states of the boson representations of 5p(12, R): 

E((N, T); KLM; T ) = aN + bN 2 + a 3 T(T + 1) + (3 3 L(L + 1) + a x T 2 (12) 

The energy of the ground state (jH)) of the system is obviously 0. 

3 Application of IVBM for the description of 
the ground state and octupole bands ener- 
gies 

In this paper we modify the earlier application of the IVBM for the 
description of the first excited even and odd parity bands in order to reach 
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much higher angular momentum states in both band types. We will apply the 
model to even- even deformed nuclei, which exhibit a low-lying negative par- 
ity band next to the ground band traditionally considered to be an octupole 
band pQ. In order to do this we first have to identify these experimentally 
observed bands with the sequences of basis states for the even representation 
of Sp(12,R) given in Table 1. We choose the SU(3) multiplet (0,//) for a 
description of the ground band, whereas for the octupole band the SU(3) 
multiplet (2,/i — 1) is used. In terms of (N,T) this choice corresponds to 
(N — 2fi,T — 0) for the positive (K n = 0+) and (N = 2fi + 2, T = 1) for the 
negative [K 1 " = CT) parity band, respectively. 

3.1 Yrast bands 

In this way, in the framework of the symplectic extension of boson represen- 
tations of number preserving U(6) symmetry we are able to consider all even 
eigenvalues of the number of vector bosons N with the corresponding set of 
pseudospins T. 

This approach is based on the fact that the energies (|12|) increase with 
increasing of N. We define the energies of each state with given L as yrast 
energy with respect to N in the two considered bands. Hence their minimum 
values are obtained at N = 2L for the ground band, and N = 2L + 2 for 
the octupole band, respectively. So for the description of the ground band 
our choice corresponds to the sequence of states with different numbers of 
bosons, iV = 0, 4, 8, ... . and pseudospin T = in the column labelled T = 
of Table 1. Respectively for a description of the negative parity band, we 
choose the set of states with quantum numbers N — 8, 12, ... and T = 1 from 
the same column T = . Since these quantum numbers uniquely define the 
SU(3) multiplets, which reduce to the corresponding values of the angular 
momenta L, the ground band belongs to the SU(3) multiplet (0, y) and the 
octupole band to (2, y — 1). In the so defined SU (3) representations for each 
N the maximal values of L appear for the first time (see Table 1.). 

According to the correspondence identified above between the basis states 
T = and the experimental data on the ground and octupole bands, the 
last term in the energy formula (fT2*Jl vanishes. The phenomenological model 
parameters a, b, a^, and (3 Z are evaluated by a fit to the experimental data. 
Their values obtained for some even-even deformed nuclei belonging to light 
actinides and rare earth region are given in Table 2. The second column gives 
the numbers of the experimental states used in the fitting procedure. 



8 



The comparison between the experimental spectra and our calculations 
using the values of the model parameters given in Table 2. for the ground 
and octupole bands of the nuclei Ra 22A , Th 22e , Sm 152 and Yb 168 is illustrated 
in Figure 1. All experimental data are taken from [To] . 

The agreement between the theoretical values obtained with only four 
model parameters and the experimental data for all the nuclei under 
consideration is rather good. 

Applying the yrast conditions relating iV and L the energies (|T2~j) for 
two considered bands can be rewritten as: 

E(L) = (3L(L+1) + ( 1 + V )L + C (13) 
The new free parameters j3, 7, 77, and £ are related to the previous ones as 

follows: 

p = 4b + p 3 , 7 = 2a -46, 77 = 80, £ = 2a + 46 + 2a 3 - (14) 

The values of j3 and 7 can be determined only from a fit to the posi- 
tive band energies, while r\ and £ are estimated from the negative ones, 
respectively. The values of the parameters (|14|) determine the behavior of 
the energies of the two bands and their position with respect to each other. 
In some cases ( 232 Th 234 U 236 U 238 U) the two bands are almost parallel. The 
shift between them depends on the parameter £. When they are very close 
they interact through the L -dependent interaction with a strength 7 + rj. 

As a result of our theoretical assumptions we obtained a simple formula 
for the energy levels. From (|13|) we can see that eigenstates of the first 
positive and negative bands consists of rotational L(L + 1) and vibrational 
L modes. The rotational interaction is with equal strength (3 in both of the 
bands. The obtained values of the parameter r/ are always negative, which 
means that the negative parity band is less vibrational than the positive one. 

3.2 The staggering 

In the collective rotational spectra of deformed even-even nuclei in this mass 
region some fine structure effects as back-banding and staggering behavior 
are observed . Odd-even staggering patterns between ground and octupole 
bands have been investigated recently [7j. In order to test further our model 
we applied on the energies the staggering function defined as |16j : 
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Stg(L) = 6A£(L)-4A£(L-l)-4A£(L+l)+A£(L+2)+A£(L-2), (15) 



where AE(L) = E(L) — E(L — 1). This function is a finite difference of 
fourth order in respect to AE(L) or of fifth order in respect to energy E(L) 
and is characteristic for the deviation of the rotational behavior from that 
of the rigid rotor. The calculated and experimental staggering patterns are 
illustrated in Figure 2. One can see a good agreement with experiment, 
as well as the reproduction of the "beat" patterns of the staggering behav- 
ior. They occur in the region where the interaction between the two 
considered bands is most strong or they cross. The correct reproduction 
of the experimental staggering patterns is due to the interaction term r]L 
in (|13|) between the positive and negative parity bands, which is a result of 
the introduced notion of yrast energies in the framework of the symplectic 
extension of the IVBM. 



4 Conclusions 



We have applied the Interacting Vector Boson Model for the description of 
the ground and octupole bands in some even-even rare earth and actinide 
nuclei up to very high spins. In spite of the simplicity of the model without 
introducing additional degrees of freedom we are able to describe both 
positive and negative parity bands. This is due to the specific definition 
of the states parity depending on the pseudospin quantum number T. 

The successful reproduction of the experimental energies and of their 
odd-even staggering was achieved as a result of their consideration as yrast 
energies in respect to the number of phonon excitation N that build the 
collective states. The introduction of this notion was possible, as we ex- 
tended the IVBM to its symplectic dynamical symmetry Sp(12,R), which 
allows the change of the number of bosons that are the building blocks of the 
model Hamiltonian. Nevertheless the Hamiltonian remains with only few 
phenomenological parameters and is still exactly solvable. Through the alge- 
braic properties of the dynamical symmetry chain relations between SU(3) 
and U(2) quantum numbers are established. Combining these relations with 
the notion of yrast energies the physical meaning of each term of the Hamil- 
tonian is clarified. In the rotational limit of the model in addition to the 
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rotational character of the considered bands an purely vibrational mode is 
appearing, which introduces also some interaction between them. This is the 
reason for the reproduction also of the fine effect of the structure of these 
bands. The obtained physically meaningful results are also simple and easy 
for use and they permit the application of the model to larger class of nuclei 
than the purely rotational ones. 

The symplectic extension of the Interacting Vector Boson Model permits 
a richer classification of the states than its unitary version and gives the 
possibility for a further consideration of other collective bands. In general 
the model proves appropriate for the description of diverse nuclear structure 
problems. 

5 Acknowledgements 

The authors are grateful for fruitful discussions and help on the subject of 
this paper to professors J. P. Draayer, J. Cseh and D. Bonatsos. This work 
was partially supported by Bulgarian Science Committee under contract 
number $ — 905. 

References 

[1] P. A. Butler and W. Nazarewicz, Rev. of Mod. Phys. 68 (1996) 349. 

[2] F. Iachello, A. Arima, "The Interacting Boson Model", Cambridge Uni- 
versity Press, Cambridge, 1987. 

[3] J. Engel and F. Iachello, Phys. Rev. Lett. 54 (1985), 1126; J. Engel and 
F. Iachello, Nucl. Phys. A472 (1987) 61; V. Zamfir and D. Kusnezov 
Phys. Rev. C 63 (2001) 054306. 

[4] Alonso et al, Nucl. Phys. A586 (1995) 100. 

[5] A. A. Raduta and D. Ionescu, Phys. Rev. C 67 (2003) 044312. 

[6] N. Minkov, S. Drenska, P. Raychev, R. Roussev, and D. Bonatsos, Phys. 
Rev. C 63 (2001) 044305. 



11 



N. Minkov, S. Drenska, Prog. Theor. Phys. Suppl. 146 (2002) 597. 

A. Georgieva, P. Raychev, and R. Roussev, J. Phys. G: Nucl. Phys. 8 

(1982) 1377. 

P. Raychev, and R. Roussev, J. Phys. G: Nucl. Phys. 7 (1981) 1227. 

N. Minkov, S. Drenska, P. Raychev, R. Roussev, and Dennis Bonatsos, 
Phys. Rev. C 55 (1997) 2345. 

A. Georgieva, P. Raychev, and R. Roussev, J. Phys. G: Nucl. Phys. 9 

(1983) 521. 

A. Georgieva, P. Raychev, and R. Roussev, Interacting Vector Boson 
Model and Negative Parity Bands in Actinides, Bulg. J. Phys. 12 (1985) 
2 147. 

H.Ganev, V. Garistov, A. Georgieva, Proceedings of the XXII Interna- 
tional Workshop on Nuclear Theory, Rila Mountains, Bulgaria, June 
15-21 2003. 

A. Georgieva, M. Ivanov, P. Raychev, and R. Roussev, Int. J. Theor. 
Phys. 28 (1989) 769. 

Mitsuo Sacai Atomic Data and Nuclear Data Ta- 
bles 31, 399-432 (1984); Level R etrieval Parameters 
http:/ / iaeand.iaea.or.at / nudat /levform.html 



D. Bonatsos, C. Daskaloyannis, S. Drenska, N. Kaoussos, N. Minkov, P. 
Raychev, and R. Roussev, Phys. Rev. C62 (2000) 024301. 



12 



Table 1. 



N (T)\ T 

V / \ 




±4 


±3 


±2 


±1 


o 













(o,o) 


< 








(2,0) 


(2,0) 
(0,1) 


2 








(4,0) 
(2,1) 


(4)0) 


4 1 







(4 0) 


(2, 1) 
(0,2) 


3 

6 1 





(6,0) 


(6,0) 
(4,1) 


(6,0) 
(4,1) 
(2,2) 


(6,0) 
(4, 1) 
(2,2) 
(0,3) 


4 










(8,0) 
(6,1) 
(4,2) 
(2,3) 


(8,0) 


3 

8 2 
1 





(8,0) 


(8,0) 
(6,1) 


(8,0) 

(6,1) 
(4,2) 


(6,1) 
(4,2) 
(2,3) 
(0,4) 

















Table 2. 



Nucleus 


n s 


a 


b 


a 3 


fa 


Ra 22i 


13 


0.0119 


-0.0022 


0.0789 


0.0155 


Ra™ 


18 


0.0269 


-0.0005 


0.0226 


0.0060 


TIi 222 


26 


0.0558 


0.0000 


-0.0557 


0.0030 


Th 22A 


18 


0.0242 


-0.0011 


0.0362 


0.0100 


Th 22fi 


20 


0.0194 


-0.0009 


0.0522 


0.0094 


TK 22 « 


18 


0.0092 


-0.0020 


0.1470 


0.0138 


Th 2:i2 


29 


0.0155 


-0.0021 


0.3244 


0.0128 


U 2 ' M 


19 


0.0124 


-0.0010 


0.3608 


0.0085 


JJ236 


25 


0.0154 


-0.0010 


0.2846 


0.0086 


JJ238 


27 


0.0142 


-0.0016 


0.2851 


0.0110 


Yb W » 


41 


0.0235 


-0.0056 


0.6512 


0.0295 


Sm™ 2 


15 


0.0194 


-0.0045 


0.4290 


0.0274 



13 



2,0 - 



> 
CD 



LU 0,5 - 




— Calculation 
Exsperiment 

— Calculation 
a Exsperiment 



8 10 12 14 16 



> 
CD 



LU 1 



Th 226 






Calculation 




* Exsperiment 




Calculation 




^ Exsperiment 







5 10 15 20 25 

Angular Momentum L 




Calculation 

* Exsperiment 

Calculation 

a Exsperiment 
j , i , i , i , i , i , i , i , i_ 



2 4 6 8 10 12 14 16 18 20 



16 



12 




— Calculation 
* Exsperiment 

— Calculation 



a Exsperiment 



10 20 30 40 

Angular Momentum L 




-3 1 ' ' ' ' ' 1 1 ' 1 ' 1 1 ' 1 1 1 

5 10 15 20 25 5 10 15 20 25 30 35 40 45 



Angular Momentum L Angular Momentum L 



